Origin of the long- wavelength magnetic modulation in Ca3Co206 
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The origin of the long- wavelength incommensurate magnetic structure of Ca3Co2 06 is discussed 
considering possible inter-chains super-superexchange paths. The experimental value of the prop- 
agation vector k=(0,0,A) with A > 1 can be reproduced only if one considers the next nearest 
super-superexchange interaction. A spin-dimer analysis using the Extended-Huckel Tight-Binding 
method confirms that, despite longer interatomic Co-Co distances, the latter interaction is indeed 
much stronger. The stability of the observed structure with respect to certain commensurate states 
is discussed. 

PACS numbers: 75.50.Ee, 75.30.Gw, 75.30.Fv 



The magnetic properties of Ca3Co206 have attracted 
considerable interest in the last decade, as this compound 
is considered as an exact experimental realization of an 
Ising triangular lattice [1]. Amongst the properties long 
debated were the nature of the ordered magnetic state 
and the still ambiguous origin of the steps in the mag- 
netization measurements at low temperature [H, [sl, 
reminiscent of that seen in single- molecule magnets. 

In the Ca3Co206 crystal structure only one of the two 
inequivalent Co^+ sites, in trigonal prismatic coordina- 
tion, is magnetic with a spin-state S=2 and a large orbital 
contribution. It is well established that these magnetic 
sites are strongly coupled ferromagnetically within chains 
running along the c-axis and that adjacent chains on the 
triangular lattice are coupled by weak antiferromagnetic 
(AFM) interactions in the a6-plane [5]. However, all 
studies to date have been undertaken in the framework 
of a quasi-2D lattice, i.e considering a simple ferromag- 
netic stacking of AFM triangular planes neglecting the 
true three-dimensional nature of the exchange interac- 
tions. Such approximation, justified in the light of the 
magnetic structures reported earlier [5], has been chal- 
lenged by the recent work of Agrestini et al. [6, 7]. Us- 
ing first magnetic X-ray diffration [6] and then neutron 
diffraction [7], these authors have unambiguously estab- 
lished that the magnetic structure is in fact incommen- 
surate, with a modulation along the c-axis of very long 
periodicity (~ 1000 A). Clearly, such a structure can not 
be stabilized in the quasi-2D limit since the only mag- 
netic interaction along c is ferromagnetic. In order to 
understand the origin of the proposed ground state it is 
important to reconsider all possible exchange interactions 
in the system. 

In the present communication, I analyze the energy 
of the magnetic structure considering only isotropic ex- 
change terms and two different super superexchange 
(SSE) interactions forming helical paths between adja- 
cent chains of the triangular lattice. In this framework, 
the observed magnetic modulation cannot be explained 
by considering only the nearest neighbour (noted J2) 
AFM SSE but requires the presence of a non-vanishing 



next-nearest AFM SSE terms (J3) which can uniquely 
stabilize the experimental value of the magnetic propa- 
gation vector. A spin-dimer analysis using the Extended- 
Hukel Tight-Binding (EHTB) method, shows that J3 is 
indeed the predominant inter-chain coupling term due to 
a stronger overlap in the molecular orbital involving the 
Co 3d xz and yz state. Whilst the observed structure cor- 
responds to the first ordered state, for small values of J3 
its exchange energy is unfavorable with respect to certain 
commensurate structures with equal moments, suggest- 
ing that the former is stabilized due to a lower entropy. 

The magnetic structure reported by Agrestini et al. 
[7] corresponds to a longitudinal amplitude modulation 
propagating along the c-axis of the hexagonal cell (space 
group i?3m, hexagonal setting [8]), as illustrated in Fig. 
[TJ In agreement with previous reports, the magnetic mo- 
ments are aligned along the c-axis, fixed by the strong 
single-ion axial anisotropy [3, 4]. However, what differs 
from previous studies is the amplitude modulation or 
spin-density wave (SDW), i.e the existence of regions in 
each chain with quasi- null ordered moments. It implies 
that the magnetic configuration within a triangular unit 
(3 adjacent chains) approaches Cl=(+,-,0) only at some 
particular points along the chain (see FigH]), unlike the 
proposed Partially Disordered Antiferromagnet (PDA) 
structure [S, ilO] which maintains this configuration along 
the entire chain. At other lattice points the configura- 
tion is almost C2=(+l/2, +1/2,-1/2) (Fig. [ij) and, if the 
structure is truely incommensurate, any other intermedi- 
ate situations between CI and C2 are found somewhere 
in the crystal. Such magnetic arrangement usually origi- 
nates from competing interactions along the chains in the 
presence of axial anisotropy since, for isotropic systems, 
competing interactions would simply lead a non-colinear 
state with fully ordered moments. The propagation vec- 
tor in the hexagonal setting is \<iH = 27r (0,0, A), with 
A ~ 1.01 at 18K. A varies with temperature, suggest- 
ing a true incommensurability rather than a locking at 
a particular fractional value. The fact that A > 1 is of 
particular importance, as will be discussed in the follow- 
ing sections. 
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For describing the magnetic structure, it is actu- 




FIG. 1: (Color online) Left: Sketch of the experimental mag- 
netic structure of Ca3Co2 06 showing the longitudinal ampli- 
tude modulation for three adjacent chains running along the 
c-axis (25 unit cells are displayed). Specific values of the in- 
plane magnetic configuration (1,-1,0) and (-■|,|r'|) cire also 
shown. Middle and right: Schematic drawing of inter-chain 
super-superexchange paths between cobalt site 1 (1) and site 
2 (2) for J2 (middle) and J3 (right). The ferromagnetic intra- 
chain coupling (Ji) is also shown. 

ally more convenient to work in the primitive rhombo- 
hedral setting which will be used in the rest of the 
communication. In the rhombohedral setting of space 
group i?3m, the unit-cell dimensions are a=b=c=6.274 
A and a=^=7=92.53°. There are two magnetic Co sites 
per unit-cell, Si and S2, at respective fractional coordi- 
nates (^,|,|) and (|,|,|). In this setting, the Co-chains 
are running along the [1,1,1] direction (also the direc- 
tion of axial anisotropy), and the propagation vector is 
k=27r(y ,y ,y ). The magnetic_structure derived experi- 
mentally by Agrestini et al. [7] corresponds to a mode 
belonging to a single irreducible representation, in agree- 
ment with the theory of second-order transition, noted 
Ti in Kovalev's notation [11]. For this mode, the cobalt 
magnetic moments M on Si (Mi) and S2 (M2), in a 
unit-cell translated by IlL={^x^^y^^z) with respect ot 
the zeroth-cell, are written as follows: 

Mi(Rl) = Mcos(k-RL) 

M2(Rl) = -Mcos (k • Rl + ttA) (1) 



where M is a vector pointing along [1,1,1], whose length 
is the amplitude of the SDW. The energy of this structure 
can be calculated, in the limit of isotropic exchange in- 
teractions, considering the intra-chain ferromagnetic in- 
teraction Ji, inter-chain super-super exchange (SSE) in- 
teractions and a phenomenological single ion anisotropy 
term DS^. Given the crystallographic parameters, two 
SSE interactions, mediated through Co-O-O-Co paths, 
must be taken into account. They correspond to inter- 
chain nearest (J2) and next-nearest (J3) neighbours, at 
interatomic Co-Co distances of 5.513 A and 6.274 A re- 
spectively. Each site (Si or S2) has six neighbours con- 
nected through J2 and six neighbours connected through 
J3. The list of neighbours is given in Table [D These SSE 
interactions, form helical paths between Co sites of ad- 
jacent chains within a triangular motif, as shown in Fig. 
dl clearly competing with the ferromagnetic intra-chain 
exchange, when J2 or/and J3 are antiferromagnetic. The 



SSE Interaction 


Site i 


Site j 




J2 


Si 


S2 


(-1,0,0) 




Si 


S2 


(0,-1,0) 




Si 


S2 


(0,0,-1) 




Si 


S2 


(-1,-1,0) 




Si 


S2 


(-1,0,-1) 




Si 


S2 


(0,-1,-1) 


Js 


Sl(2) 


Sl(2) 


(1,0,0) 




Sl(2) 


Sl(2) 


(-1,0,0) 




Sl(2) 


Sl(2) 


(0,1,0) 




Sl(2) 


Sl(2) 


(0,-1,0) 




Sl(2) 


Sl(2) 


(0,0,1) 




Sl(2) 


Sl(2) 


(0,0,-1) 



TABLE I: List of super-superexchange paths for Ca3Co2 06. 
Si and S2 refer to the two Co positions in the primitive rhom- 
bohedral unit-cell, at fractional coordinate of (^,|,|) and 
(|,|,|) respectively, while Rm refer to the translation of site 
j with respect to site i. 



energy of the magnetic mode found experimentally is eas- 
ily calculated by summing on all lattice cells (Rl)- The 
total energy, taking the convention J< for AFM inter- 
actions, is decomposed into a normal (Eat) and Umklapp 
(Ei7) terms: 



En 
Eu 



^M^ (^Jicos (ttA) + 3 J2 cos (^^^ - 3 J3C0S (^^^^ + D 



2 
2N 



Ji [cos (2k • Rl) cos (ttA)] + 3J2 



Rl 



V 3 

COS (2k • Rl) cos ( 



3J3 



cos 



(2k-RL)co5f 



(2) 
(3) 



D[cos (2k -Rl) 
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When the structure is incommensurate, as found exper- 
imentahy, one needs to consider only the non-vanishing 
normal term in equation O By derivating E at with re- 
spect to A, we obtain: 

Jisin (ttA) + J2sin ~ 2Jssin (^—^^ = (4) 

The parameter A varies in the first Brillouin zone, i.e. 
A <3/2. Let's consider FM Ji (>0) and AFM J2 and J3. 
By imposing nearest interchain neighbors only (J2 <0, 
J3=0), there are solutions uniquely for A < 1, since the 
terms in Ji and J2 must have opposite signs in Eq. [H 
However, if one considers only FM Ji and J3 <0, it is 
possible to stabilize solutions with A > 1, as observed 
experimentally. In particular, reproducing the experi- 
mental value of A ~=1.01 requires J3=-0.018Ji. 
Fresard et al. [12] already singled-out that J3 is a relevant 
parameter, based on a shorter 0-0 distance than that of 
J2. This argument is often valid, as shown by the work of 
Whangbo and co-workers in a variety of insulating oxides 
Based on the crystal structure derived by neutron 
diffraction at 18K, the 0-0 distance of 2.908 A for J3 
is indeed shorter than for J 2 (2.937 A ), despite a much 
longer Co-Co distance (6.274 A instead of 5.513 A). The 
dihedral angles along the J3 Co- 0-0- Co path also maxi- 
mize the orbital overlap as shown next. The SSE energy 
can be derived semi-quantitatively by a spin-dimer analy- 
sis based on the Extended-Huckel Tight-Binding method 
using double-^ slater orbitals for the O s and p states 
and Co d states [13]. Here, each dimer 0020^2" along 
J2 and J3 is considered in turn and the exchange energy 
is directly estimated from the square of the energy dif- 
ference between the bounding and antibounding levels 
[1^. In the present case, one needs to take into account 
the multi-electron configuration (d^, S=2 ground state). 
Since the point symmetry is high (3m), only the over- 
lap between non-orthogonal orbitals is relevant for the 
calculations as explained in [13]. The spin-dimer analy- 
sis has been performed with the program CAESAR 2.0 
[3], usin^ slater parametrizations for Co and O atoms 
given in [l5|. In trigomal prismatic configuration (point 
group 3m), the Co d orbitals are splitted into one singly 
degenerate level (z^) and two doubly degenerate levels: 
(x^-y^,xy) and (xz,yz). The z^ level is fully occupied and 
does not contribute. The four unpaired electrons occupy 
the doubly-degenerated levels, the latter lying higher in 
energy. The spin-dimer analysis shows that the energy 
difference between bounding and antibounding states in- 
volving the (x^-y^,xy) orbitals are not dramatically dif- 
ferent for J2 and J3 (considered alone they would lead 
to J2 = 1.6J3). The dimer levels involving the Co (xz,yz) 
states, form two non bounding molecular orbitals (no in- 
teraction within the numerical accuracy) and a pair of 
bounding and antibounding state, both in the case of 
the J2-dimer and J3-dimer, as displayed in Fig. [2l How- 
ever the energy gap separating the two latter states is 




FIG. 2: Bounding molecular orbitals (involving the Co 
dxzAyz levels) of the dimer Co2 0i2~ corresponding to super- 
superexchange paths along J 2 and J3, calculated from the Ex- 
tended Huckel Tight Binding method. The energy diagram 
shows the corresponding values of the squared energy differ- 
ence between bounding and antibounding states (see text for 
details). 



much stronger in the case of J3 (31 meV) than for J 2 
(7.3 meV), providing a coupling through these orbitals 
about 20 times stronger for J3. This result, even though 
semi-quantitative, supports entirely the conclusions con- 
veyed previously in the analysis of the exchange energy, 
and the initial assumption of Fresard et al. [l2[ highlight- 
ing the crucial role of the next nearest SSE interaction 
in this system. 

I now turn to the the stability of the SDW structure. 
One can easily shows that it corresponds to the first or- 
dered state, which can be derived as a function of the 
propagation vector k=(x,y,z) for various sets of exchange 
integrals {Jij}. Here the anisotropy, whose effect is to 
stabilize a SDW rather than a non-colinear configura- 
tion but does not affect the propagation, is not consid- 
ered. The ground state is obtained by calculating the 
eigenvector corresponding to the maximum eigenvalues 
of the Fourier transform of the exchange-integral matrix 
Cii [16]: 

(k, {Jij}) = ^ Jij (Rm) ■ exp {-ik ■ R^) (5) 

m 

The indices i, j refer to the magnetic atoms in a primitive 
cell (Si and S2). Jij(Rm) is the isotropic exchange inter- 
action between the spins of atoms i and j in units cells 
separated by the lattice vector R^, as listed in Table H 
In our case, there are only two magnetic atoms per unit- 
cell and the two by two Hermitian matrix is simply 
written: 

Ci,(k,U,})= (^^, ^) (6) 
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where: 
A 
B 



2 J3 {cos {x) + cos (y) + cos {z)) 



Ji (1 

J2 [e^- 



J{x-\-y-\-z) 



The phase diagram was generated numericahy with the 
program ENERMAG [l7|, performing a grid-search of 
k within the first Brillouin zone for various sets of ex- 
change parameters. Ji was fixed to 1, while J2 and J3 
were varied between -5 and +5. As expected from the 
form of ^ij, invariant by permutation of x,y and z, the 
only type of propagation-vector stabilized varies along 
a line with x=y=z. The phase diagram is represented 
in Fig. [3l In two large regions, the propagation vector 




FIG. 3: (Color Online) Top panel: First ordered state of 
Ca3Co2 06 for isotropic exchange interactions, considering a 
unitary ferromagnetic intrachain coupling (Ji = l) and a wide 
range of super superexchange interchain coupling parameters 
(J2, J3). The value of propagation vector component A is 
color-coded. Bottom panel: Energy difference between a 
SDW and the commensurate structure with k=( 1/2, 1/2,0) 
shown in the inset as a function of J3/J1 or A. 

k=0 is stabilized, with either EM or AFM configurations. 
In addition, there is a large region of incommensuratbil- 
ity corresponding to k=27r(y ,y ,y ). As already derived 
from analysis of the Fi magnetic mode energy, the line 
J3=0 stabilizes only A < 1 and the experimental value 
of A = 1.01 requires J3 7^0. However, one can also show 
that an ordered commensurate (CM) phase with equal 
moments, not obtained as first ordered state, has a lower 
exchange energy than the SDW structure. The structure, 
presented in the inset of Fig. [3] propagates with a sin- 
gle k=( 1/2, 1/2,0). Every magnetic site is fully ordered 



and has four out of its six first neighbors aligned antipar- 
allel and two aligned parallel. For values of IJ3/J1I < 
1.7, the exchange energy of this structure is lower than 
the SDW, as presented in Fig. [3l Moreover, the val- 
ues of A decreases on cooling according to the resonant 
X-ray study 0, i.e. the CM structure should be increas- 
ingly favourable. This indicate that the SDW structure, 
as observed, is stabilized only thanks to a smaller con- 
figuration entropy. The presence of competing terms in 
the free energy (exchange and entropy) could explain the 
crossover regime observed at the so-called freezing tem- 
perature T~18K, temperature below which the coherence 
length of the SDW is suddenly reduced [tI]. 

In conclusion, I have shown that the long-wavelength 
magnetic modulation in Ca3Co206 originates from the 
existence of interchain AFM super-superexchanges terms 
forming helical paths running between adjacent chains 
of the triangular motif. The analysis of the energy and 
first ordered ground state in the isotropic exchange limit, 
and a spin-dimer analysis, points to the next-nearest ex- 
change interaction as the critical parameter to reproduce 
the propagation found experimentally. It is argued that 
this phase is stabilized thanks to a lower configuration 
entropy than other structure with constant moments. 
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